Field induced stationary state for an accelerated tracer in a bath 
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Our interest goes to the behavior of a tracer particle, accelerated by a constant and uniform 
external field, when the energy injected by the field is redistributed through collision to a bath of 
unaccelerated particles. A non equilibrium steady state is thereby reached. Solutions of a generalized 
Boltzmann-Lorentz equation are analyzed analytically, in a versatile framework that embeds the 
majority of tracer-bath interactions discussed in the literature. These results -mostly derived for 
a one dimensional system- are successfully confronted to those of three independent numerical 
simulation methods: a direct iterative solution, Gillespie algorithm, and the Direct Simulation 
Monte Carlo technique. We work out the diffusion properties as well as the velocity tails: large 
v, and either large —v, or v in the vicinity of its lower cutoff whenever the velocity distribution 

&Q' is bounded from below. Particular emphasis is put on the cold bath limit, with scatterers at rest, 

3 which plays a special role in our model. 
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The propagation of a single test particle (or " tracer" ) with mass m t accelerated by an external field g in a stationary 
Q ■ bath of particles with mass mb was first studied by Lorentz in 1905 [l| as a kinetic approach of conduction in a metal. 
In this model, an electric field accelerates non-interacting electrons, which diffuse among a random configuration of 
static obstacles - the atoms - acting as infinitely massive scatterers (mt/mf, — > 0), with a simple hard sphere electron- 
atom collision rule. Lorentz assumed the velocity distribution of the electrons to obey a Maxwcll-Boltzmann law with 
a small correction, and was able to give an estimate of the conductivity of metals from microscopic parameters. It 
was later shown that this perturbed Maxwcll-Boltzmann distribution cannot be a stationary solution of the problem 
|2j : the constant acceleration confers an ever- increasing amount of energy to the electron, which cannot be dissipated 
Q ■ in the described collision process. Its velocity variance consequently diverges over time, while the mean velocity (and 
O ' hence the current) stays finite and close to the Lorentz result on a certain time scale studied in 0, [I| then vanishes at 
long times. However, a stationary solution may exist if the tracer is allowed to transfer some energy to the surrounding 
rf-\ ' medium, for instance when collisions are dissipative [5j or with a finite mass ratio mt/mt,. 

£> , The case rrn, = m t was therefore subsequently considered Q, as an extension of Rcsibois' field- less "self-diffusion 

0^ ■ problem" p} , wherein one studies the diffusive properties of identical gas particles at equilibrium by following a single 
unaccelerated tracer. The bulk of ulterior developments on this modified Lorentz problem focused on asymptotic 
and relaxation properties J8l4l3j . Recently, the model was used as a paradigm for the study of microscopic entropy 
production and fluctuation relations in non equilibrium systems [14 , 1 1 51 ] , and to discuss the force- velocity relation [16[ 
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^r: \ In what follows, we will mostly concentrate on the one dimensional version of the model, the several variants 

of which addressed in the literature differ first by the choice of different bath velocity statistics <&(v) (at rest with 
$(w) = 6(v) [a], dichotomous @, LLCS, UM or Gaussian |9|, LLCJ, UM) an d second by the collision kernel accounting for 
grain interactions, that can be indexed by a continuous parameter v (hard rods with v = 1 in most earlier works, but 

k> ■ also v = for Maxwellian rods [18[ or the so-called "very hard rods" having v = 2 [13[ ). Moreover, in addition to 
5_j \ mass ratio and field strength g that play a key role, collisions between the tracer and bath particles can be inelastic 
[ljj |2(j|. The ensuing parameter space is therefore large, with associated rich behavior. At the Boltzmann equation 
level, our goal is here to provide a comprehensive view of the tracer stationary velocity distribution function, that 
in all generality will be denoted /^(v, g,a\^): it depends on three parameters (i/, g and a number a that lumps 
mass ratio and dissipation in a single constant [Tlj, [2l| ) , and is furthermore a functional of bath velocity distribution 
<fr(v). The reason for considering an arbitrary bath distribution lies in the wide spectrum of non-equilibrium steady 
states that can be achieved upon driving granular gases with different energy injection mechanisms [20l . l22| ; Stretched 
exponentials $(t>) ~ exp[— l^l'Y/z] are often encountered here. A tracer in an evolving (e.g. freely cooling) bath may 
also be of interest [23[, assuming that it can reach a state /„ with no time dependence except through $(«,£). 

Before rationalizing the behavior of /„, the model will be detailed in section UU where a convenient integral form 
will be derived from the Boltzmann equation. This provides the basis for an efficient numerical resolution algorithm, 
that will be described in section Hi CI together with the more versatile Gillespie and Direct Simulation Monte Carlo 
techniques. Thereby equipped with three different numerical schemes, we will put to the test in section Hnl analytical 
results to be derived for the velocity tails and diffusive properties, for an arbitrary bath distribution. The cold bath 
setting will be put forward and solved in Section IIVI as a relevant model onto which more general situations can be 



mapped, while particular attention will be paid in Section [V] to the low velocity limit for bounded bath distributions. 
Finally, our conclusions will be presented in Section fVTl 

II. THE MODEL, ITS INTEGRAL REFORMULATION, AND NUMERICAL RESOLUTION 

A. Statement of the problem 

We consider the one dimensional motion of a single particle of mass m t , accelerated by a constant force g > 
through a bath. The bath is made up of unaccelerated particles of mass m&, having a stationary velocity distribution 
function <j>(i>) with characteristic velocity Vb- This velocity scale is used for the nondimcnsionalization of both the 
velocity of the tracer v = v/vb and the force g = gj ' pv\ where p is the linear mass density in the bath, v and g being 
the corresponding variables with physical dimensions. 

The time dependent tracer velocity distribution is indexed by a parameter v that specifies the type of tracer-bath 
interactions, and is denoted F v {v,t). The dependence on other parameters and on $ is left implicit unless necessary. 

The dynamics under study is governed by the linear Boltzmann (or Boltzmann-Lorentz) equation 

dF.M^ej^^ v) (1) 



dt * dv 



Here, J is the collision operator 
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where the post collisional velocity v[ of a binary encounter (vi, V2) can in general be written: 

v[ =vi— (vi - v 2 )(l - a) = avi - (a - l)v 2 (3) 

Such an expression unifies elastic and dissipative collision through the dimcnsionlcss parameter 

m t — arrib 



m t +m b 



(4) 



a being the restitution coefficient which describes the inelasticity of tracer-bath collisions (elastic if a = 1, and 
inelastic if a < 1 [l9(). It is remarkable that any system with dissipative collisions may therefore be mapped to an 
elastic system with a different mass ratio m t /mb, leading to the same dynamics and stationary states for the tracer 
[2ll . |24| . The "memory-less" situation a = is of particular interest due to its simplicity: each collision with a bath 
particle erases any memory of the tracer's previous state, and this case has been thoroughly studied and solved in 
different settings p, |9|, l2l| . The original Lorentz model corresponds to a = — 1, or a = —a in its inelastic variant [5|, 
while the opposite Rayleigh limit a = 1 (of an infinitely massive tracer) was considered in [10( . 

The exponent v introduced in @ encodes different scattering behaviors [20, [25| ; The three most common models 
are v = 1 corresponding to Hard Rods, v = for Maxwell particles and v = 2 for Very Hard Rods, the latter two 
approaches being often invoked to reproduce qualitatively some hard rod properties while simplifying calculations 
[26|] . Their merits are discussed in [l0( among others, but it will be shown in later sections that the value of v 
may affect significant properties of the solution, some of which may exhibit crossovers depending on the value of v. 
Negative values of v have been shown to produce some interesting phenomenology in the Lorentz model, such as 
the remarkable runaway effect [27| , however our one-dimensional model with finite mass ratio produces well-behaved 
solutions if v > —1 and we shall restrict our study to this case. 

It should also be noted that the molecular chaos assumption underlying the Boltzmann equation approach - 
motivated by the desire to derive analytical results- prevents the occurrence of single-file diffusion [2g. However, 
when possible, we will discuss systems with a higher dimensionality, for which the molecular chaos assumption is 
justified in the low density limit. 

It has been shown [l3[ that the velocity distribution quickly relaxes to a stationary solution /„ = lim^oc F Vl on 
which we shall concentrate in the remainder: it obeys 

dfJv) f +oc 
g— — = / dvidv 2 \vi - v 2 \" fv{vi)${v2) [5(v[ - v) - S(vi - v)] (5) 



which expresses the balance between the energy received from the external field g. and the energy transfered to the 
bath through collisions. Eq. ([5]) may be rewritten as 

9^-=r v (v)[H v (v)-Uv)} (6) 



where 
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r v (v) = I du\u-v\ u ${v) (7) 



The latter quantity is the velocity-dependent collision frequency in the gas, and H v (v) is given by 

r v {v)H v {v)= fdv 1 dv 2 \v 1 -v 2 \\fAvi)^(v2)6(v , 1 -v) (8) 



which is the gain term of the collision operator, i.e. the transition rate toward velocity v through the collision process, 
integrated over all possible initial velocities (dummy variable v\). Clearly, f u (v) = H v (v) is the equilibrium solution 
for vanishing acceleration, when detailed balance entails that enterin g an d leaving fluxes are equal at each point in 
phase space. In this field- free case (g = 0), two situations discussed in [21j ensure that the tracer velocity distribution 
is directly given by the bath distribution: first, for memory-less collisions (a = 0), we have 

/„(«,0,0|*)=*(t;) (9) 

for arbitrary bath statistics, as can be checked from ((5]). Second, when the bath distribution is Gaussian $(w) = 
G(v) = (27r)" 1 / 2 exp(-w 2 /2), one finds (see also [HI) 



\/^«) (10) 

In general however, the tracer distribution differs from the bath distribution even in this unaccelerated limit, except 
in the tails (see section ITO|) . 

We will later need to consider the distribution of pre-collisional velocities, i.e. velocities sampled right before each 
collision rather than uniformly over time. This distribution is usually biased toward high velocities -which have 
enhanced collision rate as long as v > 0- and it is given by [29| 

/„(«) = r -^-f v {v) (11) 

where the mean collision frequency is defined as 

f + OO 

dvr v (v)f v (v) (12) 



B. Implicit form 

It proves convenient to recast the stationary probability density function f v (v) as an implicit integral form. To this 
end, we write the velocity of the tracer at any given time as v = V\ +gt, with Vi the velocity that was acquired during 
its last encounter with a scatterer, and gt the contribution of the driving field during the timespan of ballistic flight 
t. We next consider the conditional probability that no collision occurs during t for a tracer starting with velocity u 1; 
which involves the velocity-dependent collision frequency r v (v): 

* /AT ft* \ 

P{t\vi)= ]J [l-AT7v(w + fffcAT)] A -^ o Cxp[- / dTr v ( Vl +gT)\ (13) 

fc=i r ^° \ Jo J 

We then need the probability density for the post-collisional velocity v\ . By definition it is proportional to the gain 
term r v (v\) H„(vi) of the Boltzmann operator, which covers all collision events from which the tracer may emerge 



with velocity v\ (proper normalization of the probability densities will be enforced a posteriori) . Finally, integrating 
over all possible time-spans t > 0, the stationary tracer velocity distribution can be written 

fu(v) = / dtexp(- drr v (v - gr)j r v (v - gt)H v (v - gt) (14) 

It can be checked by direct calculation that the above function obeys the Boltzmann-Lorcntz equation ([5]) . A similar 
equation may be found for the time-dependent solution, with explicit dependence in the initial condition F v [v, 0): 

F v (v,t)= I dt'cxpl -/ drruiv-gr) ) [r v (v - gt')H v (v - gt',t-t') + F„(v - gt,0)5{t-t')} (15) 
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where H v {v,t) is defined as in ([SJ substituting F u {v,t) for f v (v). This expression is useful for studying the relaxation 
process toward the stationary solution /„(u) in the limit t — > oo [13]. The stationary solution may also be interpreted 
as an eigenfunction of a linear integral operator £„, associated with the eigenvalue 1: 



/„(«) =£„[/„](«)= / duf v {u)K v {u,v) (16) 

J — OO 

where the kernel K u (u 7 v) is identified from equation (|14p: 

K u (u,v) = l +O °dteMRu(v-9t)-Ru(v)} '"",^+T ^( ^Y-a^ ) (17) 



with 



5 Jo 
This form simplifies in the case of Maxwell particles (y = 0) 



Rv(v) = - I dur v {u) (18) 



r (v) = / du$(u) = 1 and i? (v) = - (19) 
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The kernel then reduces to a Laplace transform 
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1 f ^ _ trK fv — gt — au 



K (u,v) = - / dte-'$ ^ (20) 

If furthermore a = 0, the solution itself becomes a transform of the bath distribution (as noted in JIG . Il3j in the 
special case of a Gaussian bath $(w) = G(v)): 

/•+oo 

fo(v,g,0)= dte-^iv-gt) (21) 

Jo 

For general values of i^ and a however, no obvious simplification can be found for the kernel K u {u 1 v). Expression (|16[) 
nevertheless allows for efficient numerical solving, as we now describe. 

C. Numerical resolution of the Boltzmann-Lorentz equation 

Three different techniques have been employed to obtain numerically the tracer velocity statistics. We start by 
summarizing their main features, before testing their compatibility. 

1. Numerical integration and iterative solving 

The operator C u derives from the Boltzmann-Lorentz equation, which by construction preserves the total probability 
in phase space. £„ must therefore conserve the integral of any function upon which it is applied. This implies that 
its spectrum is reduced to two cigcnspaccs, associated with eigenvalues 1 (where the sought solution lies) and (with 



non- physical eigcnfunctions having vanishing integral). The application of £„ therefore acts as a projection on the 
eigenspace associated with the eigenvalue 1. 

Numerically, discrctizing L v into a matrix L v and starting with a random vector, the iterative application of L v and 
subsequent normalization of the vector should converge toward the physical solution. Indeed all elements of L v are 
positive, therefore its largest eigenvalue is non-degenerate by the Perron-Frobcnius theorem, associated with a positive 
eigenvector. This unique eigenvector coincides with the physical solution obtained through numerical simulations or 
determined analytically whenever possible, and this technique allows for fast computation of the solution to any 
precision, especially in the large velocity tails where other methods may become inefficient. At this point, we have 
to choose between two strategies: solving the eigenproblcm numerically for the matrix L„, or alternatively applying 
the matrix iteratively until proper convergence is achieved. We have checked that both routes are equally precise and 
efficient, for a similar computational cost. 

2. Gillespie algorithm 

The second numerical technique used relies on the Gillespie algorithm, adapted to the Boltzmann-Lorentz equation 
by Talbot and Viot [301 . The salient features are as follows. We recall from (|13j) the conditional probability that no 
collision occurs during time t for a tracer starting with the post-collisional velocity v\ 

In P{t\vi) = R v {vi+gt) - R v (vi) (22) 

This relationship may be inverted (either analytically or numerically by dichotomy, the latter being necessary for 
v 7^ 0, 2 or with a non-Gaussian bath) to obtain t{P(v\ + gt\v±)), allowing us to randomly generate collisions times 
for the tracer, after having drawn P uniformly in [0,1]. 

The velocity vi of the collision partner must then be generated with the correct weight w(v'). First, the flux of 
collisions is divided in two contributions, coming from the right and left of the tracer: 

du \vi - u\"$(u) (23) 

-oo 
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r-(vi)= du \v! - u\ v <f>(u) (24) 

The probability of a collision coming from the right (rcsp. the left) is given by r+(vi)/r(vi) (resp. r-(vi)/r(vi)). 
The direction of the next collision is thus determined using another random uniform number over [0, 1]. Depending 
on which side is chosen, w{v2) = w±(v%) given by 

r±{V2) 

The velocity v' must be drawn with this weight, which is performed using an acceptance-rejection process: Vi is 
generated according to the bath distribution $(i>2) (which significantly improves the efficiency in generating high 
velocities over a uniform generation with a cut-off), and then accepted with a probability 

Pacc{V2) = ^7 s p^: 26) 

9{V2) W±(Vo) 

with ^o the velocity for which w(vq) is maximal (which can be computed analytically for v = 0, 1, 2). 

3. Direct Simulation Monte Carlo 

We also used Direct Simulation Monte Carlo (DSMC) [31( simulations as an alternative to the Gillespie algorithm, 
in order to experiment easily on various baths and collision laws, without performing the (bath-dependent) analytical 
computations which increase the efficiency of the Gillespie method. DSMC differs from the previous approach in that 
instead of generating the correct collision times directly, a collision partner is drawn at each step for the tracer, with 
a velocity distributed according to the bath probability density function $(«2). This collision event is then accepted 
with probability 

P acc = %^ (27) 

« "max 



with $7(t>2 — Vi) = \v2 — Vi\ v and J7 max the maximal value encountered for il at this point. The time counter is 
incremented by l/f2(v 2 — vi) at each step. This procedure converges in a few hundred steps per particle toward a 
satisfactory evaluation of the correct collision frequency, as the observed velocity distribution for the tracer (over 
many realizations, and/or sampled uniformly over time) is then identical to the result of the Gillespie algorithm, or 
analytical computations whenever they can be performed. 

This method is not significantly less efficient than the Gillespie algorithm, despite its lesser specificity. It can be 
extended to represent collisions in higher dimensionality -as required for the discussion in section Till B| - by drawing 
a vector a uniformly on the unit sphere and redefining 

(28) 



0(V 2 -Vi) = cr 1 [(V 2 -Vi).<7 



with a 1 the cross-section of the tracer-bath collisions. The post collisional velocity for the tracer is then given by 



vi = vi + (1 - a) [(v 2 - vi).<t]ct 



(29) 



4- Assessing the three numerical techniques 

The implicit route allows for efficient computation of the stationary tracer velocity distribution with important 
precision. Hereafter, it will be the method used whenever a specific prediction requires unusual numerical precision 
(good sampling of the tails). This will in particular be the case of Figs. [5] and [6) Among the three techniques, DSMC 
appears on the other hand to be the more versatile, since a change in the model parameter like softness exponent 
v does not lead to any complication. Finally, Gillespie and DSMC algorithms are well suited to study the transient 
regime before the steady state is reached (see e.g. [10|, HH, [l3| ) , together with diffusive properties. Fig. Q] demonstrates 
the validity and consistency of all three techniques in two settings. 



F (v,t) 




FIG. 1. Left: Time dependent velocity distribution Fo(v,t,g — l,a — 0), in a Gaussian bath $(u) = (2tt)~ 1 ' 2 exp(— v 2 /2) 
starting from a Gaussian initial condition with half variance. Gillespie simulation results (symbols) are compared to the analytic 
solution obtained in [ly] (solid lines). Right: stationary solution fi(v), with the same parameters and initial conditions, in a 
Lorentzian bath <!>(«) = 7r _1 /(l + v ). The results from numerical integration (solid line) and DSMC simulation at large times 
(symbols) are seen to be fully consistent. 



III. VELOCITY TAILS AND DIFFUSIVE BEHAVIOR 



A. Asymptotics 



As we shall disregard the time evolution of the system, the term " asymptotic behavior" refers here to the behavior 
of the stationary tracer distribution f v {v) for extreme values of the velocity. Rewriting the implicit equation (|14[) as 



m = \ 



du exp [R v (u) — R u (v) + In H u (u)] r v {u) 



(30) 



with H y and R v defined in ([5J and (fl8|) . and assuming that both $ and / decay exponentially we may use Laplace's 
method to approximate the integral over t. It is easily seen that 

\v\ u+1 
R v {v) ~ — — when v — >• oo (31) 

while we may choose a general stretched exponential expression for the tails of the bath distribution 

\n^(v) = --\if + 0(lnv) (32) 



where /i > 0, as found in various granular systems [20. I22L l32| . Power-law tails are obtained in the limit fj, — > 0, while 
jjl — > oo corresponds to a uniform distribution with bounded support [— 1, 1]. We shall restrict ourselves to v > — 1 
as R v becomes logarithmic when j/ — > — 1 and we cannot expect an exponential decay of the integrand and apply 
Laplace's method anymore. 

We may gain intuitive insight into the large v asymptotics by noting that competing effects are at work: such 
velocities may be reached either through collision with sufficiently energetic bath particles, or through acceleration 
without collision for a sufficiently long time. The relative importance of these effects is thus determined on the one 
hand by the abundance of energetic particles in the bath, characterized by the parameter //, and on the other hand by 
the velocity-dependent collision frequency characterized by v. Consequently, if v is large enough compared to fx, the 
increase in collision frequency for high velocities is so steep that ballistic flight is interrupted before the tracer can be 
significantly accelerated. The largest velocities are therefore reached through collisions with energetic scatterers, and 
the Boltzmann equation is dominated by the gain term, entailing that the tracer distribution is pushed back toward 
the bath distribution. This may be seen as thermalization with the bath tails. If on the other hand \i is sufficiently 
large compared to v, the bath tails are relatively depleted, while ballistic flight and acceleration are less impeded by 
the collision process. This corresponds to a prevalence of the loss term in the master equation: it is more probable 
for the tracer to have a large velocity before a collision (due to acceleration) than after it. 

A more refined asymptotic analysis confirms the above qualitative features. Since the argument of the exponential 
in pop decreases as u — > — oo, either the global maximum is located on the boundary at u = v, or it corresponds to 
some local extremum u = uq such that the derivative of the argument vanishes 

-r„M + ^^M=0 (33) 

g du 

If u = uq the v dependence of the integral receives no contribution from the In H v (u) term -the gain term of the 
Boltzmann equation is negated- and we have f v (v) ~ exp[— R v (v)]. This is possible only in the positive tail : when 
v — > — oo, all the local extrema eventually exit the interval [— oo, v] and the global maximum is necessarily located at 
the boundary u = v, regardless of the parameters [i and v. 

In the latter situation (for either tail), R v disappears from (|30|) -the loss term is negated- and we have 

IM ~ H v {v) = —^ / +D ° du \u v\" /„(«) * ( V -f^L) (34) 

We search for self-consistent solutions using Laplace's method, and find 

/„(«) ~ $(Av) (35) 

as in the unaccelerated problem, where the dilation coefficient A satisfies 

h-lolMM-i)^- 1 )/" a 

A = ! and A = if fi = 1 (36) 

1 — a 1 — a 

For a Gaussian bath (for which /i = 2) it becomes A 2 = (1 + a) /(I — a) which equals mt/rrib when collisions are elastic, 
as expected from thermalization. 

The other solution to (|34|) is f v (v) ~ f v {v/a), provided that <£>(v) vanishes faster than f v (y/a) so that the stationary 
point is asymptotically close to V = au. This solution is self-consistent only if a < 0, as it then relates one tail to the 
other rather than equating the same tail in two different points. 

Combining these observations, it appears that the leading behavior for v — > +oo is simply controlled by the value 
of the integrand of (|30[) at each of the two locations u = u and u = v, and switches from field-driven (due to the 
loss term) if the first dominates to thermalized if it is the second (see Table HJ. For v — > — oo, however, the default 
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TABLE I. Asymptotic behaviors for the positive tail v -^ +co: thermalized (left column of the table) or field-driven (right 
column). The negative tail is either its image with v —¥ av in the field-driven case with a < 0, or thermalized in any other case. 
This analysis can be extended to any v > — 1, including the three usual values v — 0, 1, 2 as shown on the right-hand graph. 

behavior is thermalized, unless a < and the positive tail is field-driven. In that case, collisions may revert large 
velocities, and the left tail becomes the image of the right tail with v — > av. 

The large velocity scenario sketched in TableUis as follows. If v + 1 > //, we recover the bath behavior in exp(— w M ), 
in agreement with the intuitive argument outlined above, while if v + 1 < //, the positive velocity tail exhibits a 
exp(— v v+1 ) decay, hereafter referred to as the cold bath behavior. These expectations, disregarding subleading terms 
such as polynomial prcfactors, are fully corroborated by Fig. f2J 

\iv + l = /-)-, both behaviors correspond to the same exponent, and prevalence of one process over the other depends 
on the value of the driving field g compared to A - * 1 . For a Gaussian bath with /i = 2, the borderline case corresponds 
to the hard-rod model, as depicted in the middle picture of Fig. O upon increasing the acceleration, the positive 
tail switches continuously from the field-independent expression exp [— (Ai>) 2 /2] when gX 2 < 1 to the field-dependent 
exp [— v 2 /2g\ when gX 2 > 1. The critical field, 1/A 2 , reduces to unity in the memoryless case a = 0. 

If v + 1 > /j,, the onset of the positive tail behavior described above is given by the point where 

v u+1 ~ gX^v" =► v ~ v c = (gX^)^ 1 -^' 1 (37) 

Assuming this threshold v c is larger than the typical bath velocity Vb (set to unity), there is an intermediate asymptotic 
range 1 <C v ^C v c dominated by acceleration. The bath- like tail is thus expected to appear only when v 3> v c . If 
however v + 1 < /x, both the intermediate range and the (positive) tail are similarly field-driven, and no such clear-cut 
crossover can be observed. 

Both cases are represented in Fig. [2]: the latter appears on the right-hand side, and the former on the left-hand side, 
where the crossover from the field-driven behavior to the bath-driven tail may be seen on the curve with v c = gX 2 = 0.5 
(crosses). The intermediate asymptotics may nevertheless dominate over a large range of positive velocities, as can 
be seen for a higher value of v c = 5 (filled dots) where the thermalized tail has yet to appear. 

In the strong- field limit g — > oo, the threshold v c diverges and the field-driven behavior prevails for all v 3> 1 over 
the whole parameter space of /j, and v. As can be seen from our choice of dimensionless variables in section Hi A[ this 
is equivalent to taking the limit of motionless bath particles v\, — > while keeping v constant : this limit $(u) — > S(v) 
is the so-called "cold bath" that will be discussed in section ITVl and exhibits characteristic field-driven tails. 

The opposite limit g — > is singular as it suppresses this cold bath (field-driven) behavior altogether, rendering our 
discussion of limiting cases invalid as the only possible asymptotic behavior is the thermalized expression f v (v) ~ $(Aw) 
in both tails. 



B. Diffusive properties 

Another point of interest lies in the tracer diffusion coefficient D, which was computed in [13| for a Gaussian 
bath and a = 0, using both its Grecn-Kubo expression and its definition in the hydrodynamic diffusion mode of the 
inhomogeneous Boltzmann-Lorentz equation. These two methods were shown to give identical results for v = and 



v = 2. The second method was also applied in [12( to the case v = 1 and a = with a dichotomous bath, revealing 
the complex dependence of the diffusion coefficient on the magnitude of the field, with a minimum at finite g, while 
this coefficient D was shown to be monotonically increasing for v = and decreasing for v = 2 in the aforementioned 
discussion of the Gaussian bath [131 ] . Here we consider an arbitrary bath distribution and parameter a, and observe 
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FIG. 2. Asymptotic behavior in a Gaussian bath (fi = 2). From left to right: f2(v), fi(v) and fo(v), computed from 
DSMC simulations for a = —1/3, g — 1 (crosses) and a = 1/3, g = 3 (dots). The solid line is the cold bath behavior 
/„(«) oc exp [— v\v\ t '/(v + l)ff]i while the dashed line is the "equilibrium" asymptotic f v (v) oc <l?(Ai;). For fz{v) (left), this 
asymptotic expression describes the positive tail, although the value of the threshold v c = gX^ may be too high for the 
crossover to be easily observable (as seen here with the dots which retain the cold bath behavior for the whole velocity range 
reported). Conversely, the negative tail is thermalized by default, and appears so even in the lone case a < 0, u + 1 < (i where 
it must become cold bath-like. This asymptote is eventually reached for lower velocities, as evidenced in Fig(4] 



once again that the effects of acceleration and thermalization disentangle when v — 0. We also demonstrate that 
these results are qualitatively unchanged in higher dimensional space, as reported for other observables in [ll| . 

As shown in the appendix, letting v — allows us to relate the diffusion coefficient to the variance of the tracer 
velocity through a simple combination 



D 



(v*)- 



1 -a 



(38) 



This variance may easily be determined recursively from that of the bath, by integration over the Boltzmann-Lorentz 
equation 



(v n )-gn(v n - 1 ) = \- 



+ OC 



du f (u) 4> n {u) 



(39) 



with 4> n {u) 



1 



-f 

1 J - c 



dvv n <& 



a - 1 " _k ~ \ a — 1 

velocity distribution is centered and has finite variance 



the nth moment of the displaced bath distribution. Therefore, if the bath 



D 



TT^(^ + (T^) 



(40) 



wher 



sre (f 2 ), is the mean velocity squared in the bath, chosen equal to 1 in our dimcnsionless variables. This result is 
in agreement with [l3j in the limit a — > 0. It may be generalized to d-dimcnsional space (see Appendix A) 



D d 



A d (l + a) 



2j2 



1 



9 2 d 



A d {l 



(41) 



where A d = 2iT d / 2 /T(d/2) is the surface area of a d-dimensional hypersphere. This prediction is in very good 
agreement with the DSMC simulation data shown in Fig. [3) In order to study the diffusive behavior of the tracer 
in these otherwise homogeneous simulations, its spatial coordinate is computed as the time-integral of the velocity. 
Between two collisions, the "space counter" in the direction of the acceleration is thus incremented by vt + \gt 2 with 
v the velocity of the tracer during this free flight and t the previously determined collision time. 
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FIG. 3. The diffusion coefficient Dd for d— 1,2,3 and g — 1 for Maxwell particles in a Gaussian bath, against the parameter 
a. The symbols represent DSMC simulations, while the solid lines are obtained using equation (|41[) . 



IV. THE COLD BATH 



In our previous discussion of the velocity tails for the tracer in a general bath, we have observed the prevalence of a 
"cold bath" asymptotic behavior for certain values of the parameters, see Table Q] The cold bath refers to a particular 
limit of the general Lorentz problem, where the scatterers are motionless and their velocity distribution is represented 
by the Dirac delta function: it may be understood as the limit Vb — > as discussed in IIII Al 

This setting has been studied for hard rods v = 1 with a = [6| and in some limiting cases with a < [lll | . Inserting 
<f>(f) = 5(v) in ([7]) gives r v (v) — \v\ v and the Boltzmann-Lorentz equation takes the form of a functional-differential 
equation 



gfl{v,g,a) = \v\ l 



1 



-jTTfv[-,g,a) -f v (v,g,a) 



(42) 



which is reminiscent of relations found in other stochastic processes such as random walks and growth models on 
trees j33[. In the field- free case, this equation reverts to a class of random collision problems systematically studied 
by ben-Avraham et al. |34| . 

This system plays a prototypical role in the study of the general Lorentz problem for two reasons. First, it may 
be seen as the strong field limit (g — > oo) of any bath distribution as discussed in section IIII Al Second, there are 
two exact mappings allowing to reduce a more general problem to the exactly solvable case of a Maxwcllian tracer 
in a cold bath, with both v = and $(u) = S(v), if only either one of these conditions is verified. Indeed, given a 
cold bath distribution, we may derive the solution for any value of the collision kernel exponent from the Maxwellian 
solution v = 0. This may be seen if we define 



9h'o( v ) 



f u (v) = h 



ho 



|f+i 



- h (v) 



(43) 



(44) 



a this equation is made to coincide with (|42[) . which is associated to Maxwell 



which verifies 



With the simple rescaling a\a\ v 
particles. 

On the other hand, for Maxwell particles iv = 0) a theorem proved by Wannier 35( for the three-dimensional case 
and Eder and Posch [l0( in one dimension gives the complete (time-dependent) solution for any bath distribution $ 
as a convolution over the cold bath solution 

/+oo 
du$(u)Fo(v-u,t,g,a\5) (45) 

-oo 

This relation holds as well for the stationary state fo(v) which is under scrutiny here. 

However Wannier's convolution theorem does not hold for arbitrary values of v. which prevents from extending this 
mapping to settings where both v and $ are chosen different from the prototypical case [y = 0, $(ij) = S(v)). Yet, 
this model still preserves many significant features of the general solution, as discussed previously. 
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A. Explicit solution 



The cold bath equation (|42[) may be rewritten as the implicit equation 

-i f+oo r 

fo( v ) — — / du fo(u) exp <d(v — au) 

gJ-oo L 5 . 



(46) 



where denotes the Heavisidc function. This expression allows us to derive straightforwardly the solution for the 
" memory- less" case a — 



fo(v,g,0) = -Q(y) exp 
9 



(47) 



where it appears that the tracer velocity distribution has support on v > 0. The solution for any v ^ may be 
recovered from the aforementioned mapping, in order to study the v dependence of quantities of interest, such as the 
mean velocity (note the choice of index v — 1 to simplify this expression) 



(v) v -i(g,a = 0) 



r ( 2 -^\ (^ 



2^ 



2v 



(48) 



where T is the Euler gamma function. Thus (v) v oc g l /( v+1 > in agreement with the dimensional argument presented 
in [6j. It should be noted that, unless v = 0, this implies a breakdown of linear response: for v = 1, (v) oc g 1 ' 2 as 
discussed in [ill ]. 

We now derive the explicit solution for a ^ 0, considering Maxwellian particles (without loss of generality, as argued 
above) which verify 



9fo(v) = ±rf (^)-fo(v) 

\a\ \a/ 

The general expression of the solution may be found by taking the Fourier transform of the function 



Ms) = 



+oo 



fo(v)e* sv dv 



(49) 



(50) 



which in turn verifies 



igsfo(s) = fo(sa) -/o(s) 



(51) 



/o(s) 



fo(sa) 
1 + igs 



(52) 



Then by recurrence we may relate /o(s) to /o( sa " +1 )- As \a\ < 1, sa n+l — > when n — > oo. Furthermore /o(0) = 1 
due to the normalization of fo(v), hence 



/o (s) = lim 



fo(sa 



n+l\ 



Uk=o (! + *ff safc ) 



fc=0 



igsa 



(53) 



Taking the inverse transform (with some care as the position of the poles and the integration contour depend on the 
sign of both a and v) 



°° p -v/ga k 

/o(f) = E — — ©(«M fe ) (i - ^T 1 



fe=0 



5 a 



(54) 



3>0,j^k 



From this expression, we may identify the asymptotic behavior for \v\ — > oo as the first non-zero term in the sum 
depending on the sign of v and a. We therefore obtain the characteristic "cold bath" positive tail behavior presented 
in Table [T] when v — > +oo, 



fo{v) ex exp 



(55) 
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The negative tail behavior however differs from the situations discussed in section Hill where it was determined by 
the bath distribution. It is readily seen from (fMf that negative velocities cannot be reached if a > : in that case, 
equation ([3|) entails that the tracer velocity is reduced but never reversed by collisions with bath particles. Once it 
has become positive due to acceleration, it may not become negative again, which ensures that the stationary solution 
has support on R + only. If however a < 0, the negative tail is given by 



fo(v —t —oo,g, a < 0) ex exp 

These asymptotic expressions for |w| — > oo with a < are corroborated by numerical results in Fig. [4] 
The moment hierarchy is obtained by a straightforward integration, and its general term is found 

(v n )=nlg n b n 

where b n = l/(a,a)„ using the <j-Pochhammer notation [36[ 



(56) 



(57) 



(a,q), 



nrioa 



n = 
aq m ) n>\ 



This symbol appears in combinatorics as a generalization of the Pochhammcr symbol or rising factorial: 



lim 



x(x + l)...(a; + n — 1) 



T(x) 



(58) 



(59) 



The expression (a, a) n converges in the limit n — > oo to a standard function of a known as the Euler function (not to 
be confused with the better known Euler totient function in number theory). 

We may rewrite the cold bath solution (|54p . using the b n notation for the product over j, and recognizing the 
exponential term as the solution for a = in (|4"?| where the acceleration has undergone the rescaling g — > ga k . We 
obtain the concise formula 



f (v, g, a) = ^ Y,^ a k{k+1)/2 b k f (v, ga\ 0) 



(60) 



fe=0 



This structure is interesting in several respects: first, the sum converges very rapidly for the needs of numerical 
computation, and is found in excellent agreement with our simulation results. Second, due to the simple form of 
/o(w,<7,0), using Wannier's theorem to derive the full solution for Maxwell particles in any bath will only involve 
taking Laplace transforms of the bath distribution. An application to the Gaussian bath is given below. Finally, even 
though this structure will not be preserved for systems with arbitrary bath distribution and v ^ 0, it suggests that 
some significant traits of the general solution may be contained in the limiting case a = 0, which is more readily solved 
in any model, see e.g. [3, Hj]. However, both analytical expressions (|54|) and (|60| shed little light on the behavior of 
the cold bath solution for positive a in the limit v — > 0, which will warrant a separate investigation in section IV1 



B. Application to the Gaussian bath 



We now revisit some properties of one of the most extensively studied settings, the Maxwellian tracer in a Gaussian 
bath {y = 0, \x = 2) 0, [ljj, [l3|. This will serve to demonstrate the application of results from the cold bath limit to 
arbitrary bath distributions, which we recall is exact for v = 0. 

From Wannier's convolution theorem we deduce that the solution fo(v,g,0) in this setting is simply given by the 
Laplace transform of the Gaussian G(v), as we previously saw in (|2T|) (see also [1J]) 



fo(v,g,Q\G) = — exp 
25 



1 

V 



1 + erf 



v — g 



V2 



(61) 



Inserting this expression in the series (|60|) allows for accurate numerical computation as well as asymptotic analysis 
for a ^ 0. 

We thus recover a result that was obtained in Ref. jlOJl , where use was made of a method specific to this choice 
of bath and collision kernel: using Mehler's formula [331, ^ is possible to compute the analytical solution of the 
Boltzmann-Lorentz equation as an expansion in terms of the Hcrmitc Polynomials H n (v) 



f (v\G) = G(Xv) J2 



hi. 

V2 



bn H,„ 



Xu_ 

V2 



(62) 
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FIG. 4. Asymptotic behavior for a Maxwellian tracer in a cold bath and a Gaussian bath (a — —1/3, g = 3). The symbols are 
Gillespie simulation results, the solid line is the exact cold bath solution (|54[) . and the dashed lines represent the asymptotic 
slopes exp(— v/g) (right, see Eq. I|55|)l and exp(— v/ag) (left, see Eq. (|56[) 1. Here, as a < and the positive tail is field-driven, 
the negative tail for the Gaussian bath is also similar to the cold bath solution, although that solution would vanish for negative 
velocities if a > and the left tail would then be thermalized (see Fig. [5J . 



However, this expansion diverges for any value of the parameters, even at vanishing velocity as H n (0) = (— 1)"' 2 nl/(n/2)\ 
for even n, while H„(0) = for odd n. Consequently, the expansion must be put under a different form to allow for 
numerical computation. It is suggested in [10( that, owing to a theorem by Euler 



b n = b 00 y £ l (-l) k a nk a k W 2 b k 

k=0 

from which, interchanging the sums, one finds as expected 

oo 

fo(v,g,a) = 5oo^(-l) fe a^ k+l ^ 2 b k f (v,ga k ,0) 



(63) 



(64) 



k=Q 



Finally, we may be interested in seeing how exactly the hierarchy of moments with a Gaussian bath relates to its 
cold bath equivalent. We compute the moment-generatrix, using the connection between the Hcrmite polynomials 
and the derivatives of the Gaussian 



M(t) 



+oo 



dve tv fo(v) = exp 



F 



2A 2 



J2( k )c^Ht)(iV2xy j H / 

3=0 KJ/ 



V2\ 



with 



(65) 



c(*) = X> n ff n& " 



(66) 



from which we derive a general expression for the fc-th moment of fo(v): 



d k M 
dt k 

Lfe/2j 



(* = 0) 



fc! 



j'=o J ' 



(67) 



The term j = (which is the leading order for large g) is recognized as the corresponding moment in the cold bath 
(|5"T)) . while the following terms appear to be specific to the Gaussian setting. 
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NON-NEGATIVE STATIONARY DISTRIBUTION 



The solution for a > exhibits some peculiar properties in the cold bath configuration: the tracer velocity can be 
reduced but never reversed by a collision. If the tracer reaches a positive velocity at a given time due to the positive 
acceleration, it will subsequently never be able to reverse its direction of motion. As mentioned in section IIV A[ this 
means that the stationary distribution has support on the positive semi-axis only. This gives a new limit to consider, 
that of small velocities v — > 0, which has no equivalent in any of the settings previously considered, and may appear 
only for a bath distribution with support on an interval bounded from below. 

A. Heuristic argument 

We shall now show that the asymptotic behavior of fo(v) for vanishing velocities is of log- normal form. We start 
with a heuristic argument. In order to reach a very low velocity v, numerous collisions must happen in a very short 
time span, before the acceleration can restore the velocity to its typical scale. Each collision multiplies the velocity 
of the tracer by a, therefore a sequence of repeated collisions may be viewed as a multiplicative random walk, or 
an additive random walk on lnw: lnijj +1 ss lnwj + In a where Vi is the velocity of the tracer following collision i. 
After a random but large number of steps (collisions occurring almost instantaneously) , we expect In v to be normally 
distributed. 

This intuition may be put on more solid grounds. In any sequence that allows the tracer to reach a very low velocity, 
collision i + 1 must take place in a time interval Ati short enough that gAU <C Vi—i, lest the tracer be reaccelerated to 
its former velocity. This interval, the upper bound for the time span between two successive collisions, becomes smaller 
as lower velocities are involved, therefore we may choose a bound of the form gAti = evi with a constant factor t«l. 
At least one collision must occur during each interval At,, which has probability 1 — cxp(— wAti) ~ wAij. Starting 
from a characteristic velocity v, wc need n such collisions to reach the velocity v n — a n v, and we may associate to 
this sequence the following probability density 



Writing x = 1/ In a, we have 



/ K)oc]>Ai ? « p) 



f (v) -exp||ln 2 (w)| 



(68) 



(69) 



We cannot expect to obtain more than the leading order in \n 2 (v) from such simple considerations, but Fig. [5] shows 
that the above expansion indeed provides the dominant behavior of the stationary velocity distribution as v — > 0, 
although increasingly difficult to evidence as a increases. 
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FIG. 5. Plot of ln/o(u)/ln 2 (u) following from (|54|l . as a function of ln 2 (i)) (increasing as v - 
g = 1. To each set is associated a horizontal line marking the expected limit x/2 = (2 In a) 
may see the approach toward the limit becoming extremely slow as a increases. 



0), for different values of a and 
, see Eq. (|69p . from which we 
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B. Log-normal characteristic approximation 

We now wish to improve upon the heuristic asymptotics (|69p , through a more formal derivation. It proves convenient 
to redefine /o(s) as the Laplace transform 



+00 



Ms) = / fo(v)e- sv dv 



(70) 



which, for imaginary arguments, coincides with the previous definition of M s ) given in section |IV Al and verifies 
again 



Ms) 



M sa ) 

1 + gs 



Iterating this relation, we obtain 



/o ( - ) = /( 



ga 



n+l 



n+1 
■>+l)(n+2)/2 TT A 

fc=l 



(71) 



(72) 



The product over /c may be rewritten as (—a, a)„+i using the g-Pochhammer symbol defined in (|58j) . and it is bounded 
from above by the finite asymptotic value (—a, a)^ given as a function of a in |38f . We may thus take the limit n — > 00, 
for which l/ga n+1 — > 00. However fo(l/g) is finite whereas 



s— >+oo 



lim s/ (s) = /o(0) = 

For /o(l/g) to be finite as n -> 00, we must have 



/o 



1 



ga 



n+l 



(n+l)(n+2)/2 



Therefore, taking s = a ("+ 3 / 2 ) 5 we find 



lim s fo(sy/a/g) oc exp 



X ln 2 (s) 



with x = (In a) x < 



(73) 



(74) 



(75) 



We emphasize here that in deriving the above asymptotic expansion, everything amounts to neglecting the loss 
term in the Boltzmann equation (|49[) . Indeed, if we expect foiv) to decrease steeply for vanishing velocities, it is 
much less probable to hold velocity o< 1 than v/a, allowing us to approximate (|49[) by 



hence for the Laplace transform 



,m-\h(l) 



gsM s ) ~ M as ) 



(76) 



(77) 



which may be iterated as in (|72|l. leaving out the factor (— a,a) n +i but giving the same asymptotic condition ([74 
This equation is verified exactly by the log- normal law ^: 



gs^(s) = ^(as) 



(78) 



*(*) 



1 



: exp 



(Ins — m) 



21 



2a 2 



s\/2tt(t 
where m and a are easily determined from the previously known parameters 

a 2 = — In a > (as < a < 1 in this section) 



(79) 



(80) 
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• — In a — In q 
2 



(81) 



As discussed in Leipnik [39| and Lopez |4(j, equation ([75)) and its counterpart in velocity space ([76]) are verified 
exactly by multiple families of functions. We note in passing that the existence of these different families of solutions 
is linked with well-known indeterminacy problems surrounding the log-normal: most notable is the non-uniqueness 
of its moments, as identical sets of moments may be found in several families of distributions. This problem has 
been studied since Stieltjes' memoir "Recherches sur les fractions continues" [41|, which elucidated the fact that the 
knowledge of the whole set of moments does not determine univocally a probability distribution. 

Assuming that /o(s) is indeed of log- normal form in the limit s — > oo, we may derive the correspondin g asymp totic 
expression of fo(v) for vanishing velocities, which later receives numerical confirmation. In the spirit of |39l. |42| . this 
expression is found from an integral expression verifying (j76[) 



pk+ioc 

ip s (v) = / exp 

J k—ioo 



z(lliV 



S(z)T(z)dz 



(82) 



where S(z) is any anti-periodic function of anti-period 1: S(z — 1) = ~S(z) such that T(z)S(z) is analytic for 
z G [k, k + 1]. Leipnik shows in [39| that the choice S(z) = (27r) _1 sin(7r,z) - that we also adopt below - equates <ps(v) 
with the characteristic function (i.e. Fourier transform) of the log-normal law. Applying Laplace's method to the 
integral, we expect for large \]nv\ 



<ps(v) « exp 



a 2 z? 



zq (hi V + m) 



S(z )T(z ) 



with zq defined as the point where the logarithmic derivative of the integrand vanishes 

S'(z Q ) , r'(z ) 



= <J"zo — (In v + m) + 



S(z ) T(z ) 



(83) 



(84) 



we may therefore take z$ « (lnw + m)/a 2 as the other terms grow at most logarithmically. Letting y — ]n.v 4- m 
suggests the asymptotic behavior 



/o(v = e^ m )ocexp 



y 

2a 2 



S(y/a 2 )T(y/a 2 ) 



(85) 



fo(v) ex exp XVy-z + 1 - m (xy)J for y = l n («/Va.9) 



-oo 



(86) 



where we used x^ nv ^ 1 an d Leipnik's choice S^ 1 (z) = 2r(z)T(l — z) j39|. We thus recover the leading order of 
the heuristic expression (|69|) . with additionally a sublcading contribution (^lnu) - ^ in the probability distribution, 
stemming from the gamma function. Although subleading, this correction is nevertheless important: the complete 
structure of expression (|86p is required so that, upon insertion in the Boltzmann equation without loss term, both 
sides may be matched for vanishing v. Indeed, it verifies relation (|76[) to leading order in v/hxv whereas our first 
log-normal approximation (j69[) does not. 

This new asymptotic expression is clearly seen in Fig. [6] (left) to approximate the solution fo(v) much more closely 
and on a wider range of velocities than (|69)) . Furthermore the right-hand graph demonstrates the same behavior 
for a different non-negative bath distribution: a half-Gaussian, which strictly vanishes for v < 0. This tail behavior 
appears for a > as long as there is a lower bound to the bath particle velocity. In addition, any distribution 
<&(i>) ^= 5(v) but bounded from below with support in [0, oo[ leads to a decrease of tracer velocity probabilities (see the 
half-Gaussian results in Fig. [51 well below their cold bath counterpart), shifting the tail by a multiplicative constant 
without affecting its functional form. However, the range of validity of this asymptotic expression is relegated further 
in the tail when a increases. 



VI. CONCLUSIONS 



Within a generalized Boltzmann-Lorentz framework, we have studied the stationary state of a tracer particle 
immersed in a bath with arbitrary velocity statistics $(i>). The tracer (say a charged particle) is accelerated by 
an external field, that does not act on the bath (say made up of neutral particles). In addition to <J>, the model 



17 




10 



10' 



10 



10 



10 



10 



10 



in 



15 



20 



25 



■:{{") 



Cold bath 




FIG. 6. Left: Validity of the asymptotic expressions. The exact cold bath solution fo(v) [solid lines, corresponding to (|54|) ] 
against the predictions (|69|) (dashed lines) and (|86[) (dots). Here, a — 0.01,0.1 and 0.5 from left to right, and g — 1. As a 
increases, the predicted asymptotic behavior appears further in the tail, and the expansion (|86|l becomes significantly more 
precise than our first heuristic guess (|69[1 , Right: Close-up on the data for a = 0.01 (solid line) and comparison with the solution 
corresponding to a half-Gaussian bath $(u) = 0(«)(27r) -1 ' exp(— v /2) (dashed line), other parameters being identical. The 
same asymptotic form (|86[) (dots) holds in both cases. Similar behavior is expected for any bath distribution with a lower 
velocity bound, provided g > and a > 0. 



is specified by the field intensity g, a parameter a accounting for material properties that combines mass ratio 
(tracer over bath particle) together with collisional dissipation, and an index v quantifying the softness of particle 
interactions. Maxwellian, hard, and very hard rods are thereby embeded in a unifying approach, that lends itself 
to analytical progress and numerical investigation. Three independent simulation techniques were used to solve the 
Boltzmann-Lorentz equation. We have performed an asymptotic analysis of the high-energy tails. We have evidenced 
a winner-takes-all competition between the two processes that leads to a non equilibrium stationary state: on the one 
hand the acceleration during the variable timespan of ballistic progress, which feeds energy into the system, and on 
the other hand, the collision with scatterers that lead to dissipation. Our analysis has revealed that, depending on 
the parameter and velocity range that we consider, either of these effects dominates and completely determines the 
local shape of the distribution. Therefore, significant portions of this distribution are almost independent either from 
the bath or from the acceleration and collision kernel. The latter effect may provide a way to probe non-gaussianitics 
in the bath using a more visible tracer and charting out its large velocities behavior; Such a behavior requires that 
the scattering exponent v is large enough once 3> has been chosen, or conversely that the bath distribution is not too 
cold -i.e. not too peaked around the origin- once the scattering exponent is fixed. Whereas most of the analysis was 
performed for a one dimensional system, we have, for diffusive properties, also considered driven discs or spheres in 
higher dimension. 

Motivated by the existence of bath-independent properties, particular emphasis was put on an already introduced 
setting, but previously unsolved in its general formulation: the "cold gas" with vanishing temperature [i.e. $(u) = 
5(f)] ■ Its interest as an approximation for more realistic models was addressed. We have in particular shown that 
even when the cold bath asymptotics does not prevail, an intermediate asymptotic cold bath regime may exist, when 
the tracer velocity is significantly larger than the bath characteristic velocity but lower than the model-dependent 
threshold v c . This setting also allowed us to exhibit an interesting asymptotic behavior at low velocities for a tracer 
in a bath with a non-negative velocity distribution -e.g. a heavy intruder falling in a static bath or along a stationary 
stream of particles. 
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APPENDIX: COMPUTATION OF THE DIFFUSION COEFFICIENT 

Let us define the velocity autocorrelation function T(t) as 

r(t) = ((«(*)-(«» («(o) -(«») (87) 

where ( . ) denotes the mean over the stationary distribution f v {v) and v(t) is the velocity of a given realization of the 
tracer at time t. A standard relation p, [43| connects T(t) to the diffusion coefficient 

1 d 



^^Jti^-^ 2 ) 



lim((v(t)-(v}) dt'(v(t')-(v}) 
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dtT(t) 
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where use was made of the fact that both v(0) and i>(oo) are sampled according to f v (v), and it was assumed that 
the above limits exist. In a derivation similar to the one used in [431 ] . this function T(i) is expressed as 

T(t) =fdv f v (v )fdvF v (v,t\v ) (v ~ (v)) (v - (v)) (89) 

where -F y (u,i|wo) is the conditional velocity distribution of the tracer at time t knowing it had velocity vq at time 0, 
i.e. it is the time-dependent distribution F v {v,t) with the initial condition F v (v,0) = S(v — vq), which may also be 
represented using the implicit formulation (|15[) . Wc may thus define the auxiliary function 



N(v,t)= I dv f„(vo) F v (v,t\v ) Oo - (v)) (90) 

which fulfills the initial condition 

N(v,0) = (v-(v))U(v) (91) 



r(£) = / dv(v-(v))N(v,t) (92) 



and 



Due to the linearity of the Boltzmann-Lorentz equation, N(v, t) follows the same equation as F v (v, t) as can be seen 
from Eq. (|90|) . Therefore, the determination of T(t) using ([92]) amounts to computing the integral and first moment 
of the solution of the Boltzmann-Lorentz equation with the non-physical initial condition (|iTl"T) . As this equation 
conserves the normalization, 

J dv N(v, t)= I dv N{v, 0) = (93) 

thus 

T(t) = J dvvN(v,t) (94) 



Finally, if v — 0, this time-dependent first moment can be computed directly from the Boltzmann equation 

Op p Q TIT p p 

— = -g / dvv— h / dvidv 2 N(vi)®(v 2 ) dvv [S(v[ - v) - 8(v\ - v)] (95) 

The first term on the right-hand side vanishes, while v[ — v\ = (a — 1) (v\ — v 2 ) and therefore 

^ = -(l-o)r(t) (96) 

r(t) = r(o) e -( 1 - a )* (97) 

from which the diffusion coefficient follows 

a.roEUW-W' ,98) 

1 — a 1 — a 

We have thus related D to the variance of the tracer stationary velocity distribution in the case of Maxwell particles 
with an arbitrary bath. 

Furthermore, the approach can be generalized to higher dimensions, where the collision law reads 

vi=vi + (l-a)[(v a -vi).*]£ (99) 

so that equation (|95"|) becomes 



? = - /dvv(g.V v ) N(v,i)+ / > dv 1 dv 2 d ( rN(v 1 )$(v 2 ) f dvv[6{v' 1 -v)-6(vi-v)} (100) 
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Then, if we define ei2 the unit vector along (v 2 — Vi), 9 the angle between a and ei2, and Cl the solid angle 

dfr[(v 2 - V!).ct]o- - / dh\ v 2 - vi | cos 2 (<?) e 12 = C(d) (v 2 - Vl ) (101) 



C(d) = / dft[l-sin 2 (6>)] = A d 



J*d9 sm d 
1 ~ jjf d9 sin d - 2 



d 



(102) 



where A^ is the area of the d-dimensional unit hypcrsphere. We may then compute the first and second moments: 
letting g = ge x 



v = (v x ) e x = 



gd 



A d {l-a)) 



(103) 
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X 7 l + a X /b 1 



2ffd (f g ) 
-a 2 A d 



hence 



r(o) = (v 2 ) - (vy 



1 -a 



< A + 



9 2 d 2 



l + a x /b A 2 (l-a 2 ) 
which gives the following relation for the diffusion coefficient in d dimensions in a bath with unit temperature: 



(104) 
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